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H ! Abstract 

Normal surface theory is a central tool in algorithmic three-dimensional topology, and the 
enumeration of vertex normal surfaces is the computational bottleneck in many important al- 
gorithms. However, it is not well understood how the number of such surfaces grows in relation 
to the size of the underlying triangulation. Here we address this problem in both theory and 
practice. In theory, we tighten the exponential upper bound substantially; furthermore, we 
construct pathological triangulations that prove an exponential bound to be unavoidable. In 
practice, we undertake a comprehensive analysis of millions of triangulations and find that in 
general the number of vertex normal surfaces is remarkably small, with strong evidence that 
our pathological triangulations may in fact be the worst case scenarios. This analysis is the 
first of its kind, and the striking behaviour that we observe has important implications for the 
feasibility of topological algorithms in three dimensions. 



> ■ 1 Introduction 

00 

Geometric topology is an inherently algorithmic subject, with fundamental questions such as the 
iy-} ■ homeomorphism problem (find an algorithm to determine whether two given spaces are topologically 

equivalent) and the identification problem (find an algorithm to determine the topological name 
and/or structure of a given space). Three-dimensional topology is of particular interest, since 
in lower dimensions such problems become trivial [27], and in higher dimensions they become 
unsolvable [25]. 

Throughout this paper we restrict our attention to closed 3-manifolds. In essence, a closed 
^ 3-manifold is a compact 3-dimensional topological space that locally looks like R 3 at every point. 

Much recent progress has been made on algorithms in 3-manifold topology. For example: 



o 



Rubinstein gave an algorithm in 1992 for recognising the simplest of all closed 3-manifolds, 
namely the 3-sphere [32, 33]; this algorithm has been refined several times since [8, 19, 34]. 

In 1995, Jaco and Tollefson gave an algorithm for breaking a 3-manifold down into a connected 
sum decomposition (essentially a topological "prime decomposition") [21]. 

Perelman's proof of the geometrisation conjecture in 2002 finally resolved the general home- 
omorphism problem for 3-manifolds, completing a programme initiated decades earlier by 
pioneers such as Haken [15] and Thurston [36]. The full homeomorphism algorithm is a fusion 
of diverse and complex components, including both the 3-sphere recognition and connected 
sum decomposition algorithms above. 



*A more detailed version of this paper will be available as Extreme cases in normal surface enumeration (in 
preparation) . 
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A recurring theme in these algorithms (and many others) is that they rely upon normal surface 
theory, a tool that allows us to convert difficult topology problems into simpler linear programming 
problems. In particular, we can search for an interesting surface within a 3-manifold by (i) con- 
structing a high-dimensional polytope, (ii) enumerating the "admissible" vertices of this polytope, 
and then (iii) testing each admissible vertex to see whether it encodes the interesting surface that 
we are searching for. 1 

The concept of an "interesting surface" depends on the application at hand. For instance, in the 
connected sum decomposition algorithm we search for embedded spheres within our 3-manifold; in 
other algorithms we might search for non-trivial embedded discs [14] or embedded incompressible 
surfaces [17]. However, in all of these applications the high-dimensional polytope and its admissible 
vertices remain the same. That is, the polytope vertex enumeration problem is a common component 
for all of these topological algorithms and many others besides. 

Furthermore, this common vertex enumeration problem is in fact the computational bottleneck 
for many of these algorithms [8, 10]. It is therefore important to improve the efficiency and un- 
derstand the complexity of this vertex enumeration problem, since any improvements or results 
will have a widespread impact on computational 3-manifold topology as a whole. This impact also 
extends beyond three dimensions — for instance, in J^-manifold topology, to understand whether a 
given triangulation represents a 4-manifold we require all of the complex machinery of 3-sphere 
recognition as discussed above. 

In general, polytope vertex enumeration is difficult. The general problem is known to be NP- 
hard [12, 22], and the range of available algorithms is matched by a range of pathological cases that 
exploit their weaknesses [2]. However, in our context we have two advantages: 

• We are not dealing with an arbitrary polytope, but rather one that derives from the machinery 
of normal surface theory; this polytope is known as the projective solution space. Such 
polytopes have additional constraints on their dimensions and the equalities and inequalities 
that define them. 

• We do not need to enumerate all vertices of the polytope, but only the admissible vertices. 
These are the vertices that satisfy an additional family of non-linear constraints, known as 
the quadrilateral constraints. 

These contextual advantages can be exploited in vertex enumeration algorithms with great 
success; see [7, 8, 9, 37] for details. Nevertheless, the enumeration problem remains a difficult one. 
In particular, Agol et al. [1] show that determining knot genus — yet another problem that employs 
normal surface theory — is in fact NP-complete. 

In this paper we concern ourselves with the complexity of the enumeration problem. More 
specifically, we focus on the number of admissible vertices of the projective solution space, which 
we denote by a. This quantity is important for the following reasons: 

• The admissible vertex count a gives a lower bound for the time complexity of vertex enumer- 
ation. Moreover, for the quadrilateral-to-standard conversion algorithm (a key component of 
the current state-of-the-art enumeration algorithm), there is strong evidence to suggest that 
the running time is in fact a low-degree polynomial in a [7]. 

• Each admissible vertex corresponds to a surface in our 3-manifold upon which we must run 
some subsequent test. For some problems (such as Hakenness testing [10, 17]) this test is 

1 Some other algorithms (such as knot genus [16] and Heegaard genus [24]) replace step (ii) with the more difficult 
enumeration of a Hilbert basis for a polyhedral cone, yielding what arc known as fundamental surfaces. 
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extremely expensive, and so the number of admissible vertices becomes a critical factor in the 
overall time complexity. 

The input for a typical normal surface algorithm is a 3-manifold triangulation, formed from n 
tetrahedra by joining their An faces together in pairs. We call n the size of the triangulation; not 
only does n represent the complexity of the input, but both the dimension and the number of facets 
of the projective solution space are linear in n. 

The growth of a as a function of n is currently not well understood. The only general theoretical 
bound in the literature is a < 128™, proven by Hass et al. [16]; in the special case of a one-vertex 
triangulation this has been improved to a £ 0(15 n ) [9]. Very little is known about the growth 
of a in practice, though initial observations suggest that a is in fact far smaller [7]. For example, 
in the proof that the Weber-Seifert dodecahedral space is non-Haken (one of the first significant 
computer proofs to employ normal surface theory), a "typical" triangulation of size n = 23 is found 
to generate just a = 1751 admissible vertices [10]. 

In this paper we shed more light on the growth of a, including new theoretical bounds and 
comprehensive practical experimentation. Following a brief outline of normal surface theory in 
Section 2, we present the following results: 

• In Section 3 we show that a G 0(<fi 7n ), where <\> is the golden ratio (1 + \/5)/2. This tightens 
the general theoretical bound on a from 128™ to just over 0(29"). We prove this by extending 
McMullen's upper bound theorem [31] to show that any convex polytope with k facets must 
have 0(4> k ) vertices. 

We push this bound from the other direction in Section 4 by constructing an infinite family 
of 3-manifold triangulations for which a = 17™/ 4 + n/4. This yields the first known family for 
which a is exponential in n, and disproves an earlier conjecture of the author that a 6 0(2™). 
By extending this family to all n > 5 we show that any theoretical upper bound must grow 
at least as fast as ft(17™/ 4 ) ~ ft(2.03 n ). 

• In Section 5 we build a comprehensive census of all 3-manifold triangulations of size n < 9, 
and measure a for each of the ~ 150 million triangulations that ensue. We find a remarkably 
slow growth rate — for n > 5 the worst cases are precisely the infinite family above, suggesting 
that the lower limit of f)(17™/ 4 ) ~ £1(2.03™) may in fact be tight. In the average case the 
mean a appears to grow even slower, with an apparent growth rate of less than 4> n and a final 
mean of just a ~ 78.49 for n = 9. 

This analysis is the first of its kind, primarily because the complex algorithms and software 
required for such a comprehensive study did not exist until very recently [5, 7]. Previous 
censuses have focused on restricted classes of triangulations (such as minimal triangulations 
of irreducible or hyperbolic manifolds [5, 11, 26, 30]), and previous measurements of a have 
been for isolated or ad-hoc collections of cases [7, 10, 28]. 

Throughout this paper we work with Haken's original formulation of normal surface theory 
[14, 15]. Tollefson defines an alternative formulation called quadrilateral coordinates [37], which is 
only applicable for some problems but where the polytope becomes much simpler. In quadrilateral 
coordinates an upper bound of a < 4™ can be obtained through an analysis of zero sets [9], but 
again the growth rate is found to be significantly slower in practice. We address quadrilateral 
coordinates in detail in the full version of this paper. 
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2 Preliminaries 



Throughout this paper we assume that we are working with a 3-manifold triangulation of size n. 
By this we mean a collection of n tetrahedra, some of whose 4n faces are affinely identified (or 
"glued together") in pairs so that the resulting topological space is a 3-manifold (possibly with 
boundary). If all 4n faces are identified in 2n pairs then we obtain a closed 3-manifold; otherwise 
we obtain a triangulation with boundary, and the unidentified faces become boundary faces. Unless 
otherwise specified, all triangulations in this paper are of closed 3-manifolds. 

There is no need for a 3-manifold triangulation to be rigidly embedded in some larger space — 
tetrahedra can be "bent" or "stretched" . Moreover, we allow multiple vertices of the same tetra- 
hedron to be identified as a result of our face gluings, and likewise with edges. This allows us to 
build triangulations using very few tetrahedra, which becomes useful for computation. 




Figure 1: A 3-manifold triangulation and an embedded normal surface 

To illustrate, the left-hand diagram of Figure 1 shows a triangulation of the product space 
S 2 x S" 1 using just n = 2 tetrahedra — the back two faces of each tetrahedron are identified with 
a twist, and the front two faces of the left tetrahedron are identified directly with the front two 
faces of the right tetrahedron. All eight vertices become identified together, and the 12 edges 
become identified in three distinct classes (represented in the diagram by three different types of 
arrowhead). We say that the resulting triangulation has one vertex and three edges. 

Normal surfaces were introduced by Kneser [23], and further developed by Haken [14, 15] for use 
in algorithms. A normal surface is a 2-dimensional surface embedded within a 3-manifold triangula- 
tion that meets each tetrahedron in a (possibly empty) collection of triangles and/or quadrilaterals, 
as illustrated in Figure 2. For example, a normal surface within our S 2 x S 1 triangulation is shown 
on the right-hand side of Figure 1; as a consequence of the tetrahedron gluings, the six triangles 
and quadrilaterals join together to form a 2-dimensional sphere. 




Figure 2: Normal triangles and quadrilaterals within a tetrahedron 

There are four distinct types of triangle and three distinct types of quadrilateral within each 
tetrahedron (defined by which edges of the tetrahedron they meet). The vector representation of 
a normal surface is a collection of 7n integers counting the number of pieces of each type in each 
tetrahedron; from this vector in R 7n we can completely reconstruct the original surface. We treat 
surfaces and their vectors interchangeably (so, for instance, "adding" two surfaces means adding 
their two vectors and reconstructing a new surface from the result). 

An early result of Haken is a set of necessary and sufficient conditions for a vector to represent a 
normal surface: (i) all coordinates must be non-negative; (ii) the vector must satisfy a set of linear 
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homogeneous equations (the matching equations); and (hi) there can be at most one non-zero 
quadrilateral coordinate corresponding to each tetrahedron (the quadrilateral constraints). Vectors 
that satisfy all of these conditions are called admissible. 

Jaco and Oertel [17] define the projective solution space to be the polytope in R 7n obtained as a 
cross-section of the cone defined by (i) and (ii) above. A vertex normal surface lies on an extremal 
ray of this cone and is not a multiple of some smaller surface. The vertex normal surfaces are in 
bijection with the admissible vertices of the projective solution space; we let a denote the number 
of vertex normal surfaces, and we call a the admissible vertex count. 

The enumeration of vertex normal surfaces is a critical component — and often the computational 
bottleneck — of many important topological algorithms. This is because one can often prove that, 
if an interesting surface exists (such as an incompressible surface or an essential sphere), then one 
must appear as a vertex normal surface. See Hass et al. [16] for a more detailed introduction to 
normal surface theory and its role in computational topology. 



3 Theoretical Bounds 

As noted in the introduction, the best bound known to date for the admissible vertex count is 
a < 128 n , proven by Hass et al. [16]. We begin by tightening this exponential bound as follows: 

Theorem 1. Let (ft = (l + \/5)/2. Then the admissible vertex count a is bounded above by 0((j) 7n ) ~ 
O(29.03 n ). 

We prove this through a simple extension of McMullen's upper bound theorem [31]. McMullen 
gives a tight bound on the number of vertices for a convex polytope with k facets and d dimensions; 
we extend this here to a loose bound that covers all possible dimensions. 

Lemma 2. Let Fq = 0, F\ = 1, F2 = 1, ... represent the Fibonacci sequence, where Fi + 2 = 
Fi + \ + Fj. Then for any k > 3, a convex polytope with precisely k facets has < F k+ i vertices. 

Proof. Suppose the polytope P is <i-dimensional with precisely k facets. Then McMullen's theorem 
(taken in dual form) shows that P has at most 

k ~ ^ fk - _ fk- [^\\ , (k - 



k-d ) + \ k-d ) ^_ L ^±ijJ + U-L^jJ (1) 

vertices. 2 For even d this can be rewritten as ( k ~ a ) + (^ _ jf~ 6 ) f° r suitable integers a, b, and for 

odd d it can be rewritten as 2{^ k ~ 1 ^~ a ) for a suitable integer a. 

We now claim that ( h ~ a ) < F k for any k, a with k > 1. This is easily established for k = 1, 2, and 

the full claim follows from°the inductive step { k ~ a ) = ( fe ~*~ a ) + (^7°) = (^""M**" 2 ^ 0-1 *) < 
iVi + F fe _2 = F k . 

From here our lemma is straightforward. If d is even then the number of vertices of P is at most 
{ k ~ a a ) + ( (fc ~ 2) " 6 ) - F k + F k-2 < Fk + Fk-i = Fk+i, and if d is odd then the number of vertices is 
at most 2(( fe -p~ a ) < 2F k ^ < F k + F k ^ = F k+1 . □ 

Unlike McMullen's result, Lemma 2 is not tight. Nevertheless, it gives us a very good 3 asymp- 
totic upper bound of 0{<j) k ), which is enough to prove our main theorem. 



2 This is the number of facets of the cyclic d-dimensional polytope with k vertices [13]. 

Experimentation shows that this asymptotic upper bound of (f> k ~ 1.618 fe is close to optimal. If we maximise 
equation (1) over all d for each k = 100, . . . , 200, the maximum grows at a rate of approximately 1.613 fe . 



5 



Proof of Theorem 1. The facets of the projective solution space in W. 7n are defined by the 7n in- 
equalities x\ > 0, . . . , x-jn > 0, and so there are at most 7n facets in total. Lemma 2 then shows 
that the projective solution space has at most F-j n+ \ vertices, and so a < F-j n+ \. Using the standard 
formula F k = [cf) k /V5 + \\ it follows that a G 0((f) 7n ). □ 

It is interesting to note that Theorem 1 makes no use of admissibility — this suggests that, 
although the bound of 4> 7n is a strong improvement on 128 ra , this bound is still very loose. We 
confirm this through experimentation in Section 5. Although we only consider closed 3-manifolds 
in this paper, it should be noted that Theorem 1 and its proof apply equally well to triangulations 
with boundary, and also to the ideal triangulations of Thurston [35] . 

4 Extreme Cases 

Having tightened the upper bound from above, we now turn our attention to limiting the upper 
bound from below. We do this by building pathological triangulations for which a G 6(17 ra / 4 ) ~ 
6(2. 03 n ). This growth rate shows that an exponential upper bound on a is unavoidable, and 
furthermore disproves an earlier conjecture of the author that a G 0(2 n ). 

We begin by describing 4-blocks, which are small building blocks that appear repeatedly 
throughout our triangulations. Using these building blocks, we then construct the family of patho- 
logical triangulations X±, X2, ■ ■ ■■ 

Definition (4-block). A 4-block is a triangulation with boundary, built from the four tetrahedra 
Ai, A2, A3, A4 using the following construction. 




Triangular pillow 

Figure 3: The two-tetrahedron triangular pillow at the centre of a 4-block 

We begin by folding together two faces of A l5 and then wrapping A 2 around the remaining two 
faces as illustrated in Figure 3. This forms a triangular pillow with three vertices, three boundary 
edges, two internal edges, and two boundary faces. 

Next we fold together two faces of A3 and two faces of A4, as illustrated in the leftmost column 
of Figure 4. To finish, we join the pillow to both A3 and A4 as illustrated in the central column of 
Figure 4 — the upper face A\B\A2 of the pillow is glued to the lower face A^B2A^ of A3, and the 
lower face A\B\A2 of the pillow is glued to the upper face A4B3A4 of A4. 

The final result is shown in the rightmost column of Figure 4, with three boundary vertices 
and one internal vertex. The triangular pillow is buried in the middle of this structure, wrapped 
around the internal vertex; for simplicity the two edges inside the pillow are not shown. 

Definition (Pathological triangulation X^). For each integer k > 1, the pathological triangulation 
Xk is constructed from n = Ak tetrahedra in the following manner. 

From these 4k tetrahedra we build k distinct 4-blocks, labelled B\, . . . ,Bk- Within each 4-block 
Bi we label the three boundary vertices Pi,Qi,Ri, where Pi sits between both boundary triangles 
as illustrated in Figure 5. 
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Initial pieces Gluing instructions 



Figure 4: Building a 4-block from two tetrahedra and a triangular pillow- 




Figure 5: Building the pathological triangulation X\. from k distinct 4-blocks 



For each i = 1, . . . , k we join blocks B{ and as follows (where Bk+i is taken to mean B±). 
Triangle P{P{Ri is joined to triangle Qj+iPj+iPj+i; note that this is "twisted", not a direct gluing, 
since it maps P; L Qi+i and There are in fact two ways this gluing can be performed 

(one a reflection of the other); we resolve this ambiguity by orienting each block consistently, and 
then choosing the gluing that preserves orientation. 

An effect of these gluings is to identify all of the Pi, Qi and Ri to a single vertex, so that X^ 
has k + 1 vertices in total (counting also the k internal vertices from each original block). 

It is not clear that each X^ is a 3- manifold triangulation (in particular, that X\. looks like M 3 
in the vicinity of each vertex) . The following sequence of results proves this by showing that every 
Xk is in fact a triangulation of the 3-sphere. 

Lemma 3. A 4-block is a triangulation of the 3-ball (i.e., the solid 3- dimensional ball), with a 
boundary consisting of two triangles in the formation shown in Figure 6. 

Proof. This is evident from the construction in Figure 4. It can also be verified computationally 
using the software package Regina [3], which implements 3-sphere and 3-ball recognition [4]. □ 

Lemma 4. Let 71 and T2 each be triangulations of the 3-ball with boundaries in the formation 
shown in Figure 6. If we identify one boundary triangle of 71 with one boundary triangle of 7i 
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Triangle XYX 







Triangle XZX 




Figure 6: A 3-ball whose boundary consists of two triangles 



under any of the six possible identifications, the result is always another triangulation of the 3-ball 
with boundary in the formation shown in Figure 6. 

Lemma 5. Let T be a triangulation of the 3-ball with boundary in the formation shown in Figure 6. 
If we identify the two boundary triangles under any of the three possible orientation-preserving 
identifications, the result is always a closed 3-manifold triangulation of the 3-sphere. 

Proof. Both of these results are essentially properties of 3-manifolds, not their underlying triangu- 
lations — if they hold for any selection of triangulations 71 , T2 , T then they must hold for all such 
selections. We verify these results using Regina by choosing 4-blocks for our triangulations and 
testing all six/three possible identifications. □ 

Since each X^ is built by joining together 4-blocks along boundary triangles in an orientation- 
preserving fashion, the following result follows immediately from Lemmata 3-5. 

Corollary 6. For each k > 1, X^ is a closed 3-manifold triangulation of the 3-sphere. 

We turn our attention now to counting the vertex normal surfaces for each triangulation X^. 
Recalling that k = n/4, the following result shows that for these pathological triangulations we 
have a G 9(17 n / 4 ) ~ 6(2.03"). 

Lemma 7. For each k > 1, Xy. has precisely a = 17 k + k vertex normal surfaces. 

Proof. Consider a single 4-block with boundary vertices labelled P, Q, R as before, and let S denote 
the internal vertex. Define a, f3 and 7 to be small loops on the 4-block boundary surrounding P, 
Q and R respectively, as illustrated in Figure 7. 



Using the software package Regina, we can construct the projective solution space for this 4- 
block. There are 17 admissible vertices in total, corresponding to 17 vertex normal surfaces: one 
with empty boundary, and 16 whose boundary consists of some combination of a, (3 and 7. These 
surfaces are summarised in Table 1, and we label them a, b, . . . , q as shown. 

It is important to note that a, b, . . . , q are all compatible; that is, no combination of their vectors 
can ever violate the quadrilateral constraints. 4 This is an unusual but extremely helpful state of 
affairs, since we can effectively ignore the quadrilateral constraints from here onwards. 

4 This is because, within each tetrahedron, we observe that two of the three quadrilateral types never appear 
anywhere amongst the surfaces a, b, . . . , q. 




Figure 7: The curves a, (3, 7 on the boundary of a 4-block 



8 



Label 


Boundary 


Description 


a 




— 


Small sphere around internal vertex S 


b 


a 




Small disc around boundary vertex P 


c 






Small disc around boundary vertex Q 


d 




7 


Small disc around boundary vertex R 


e 


a 




Tube from P to S, closed around S 


f 






Tube from Q to S, closed around S 


g 




7 


Tube from R to S, closed around S 


h 


a 


+ P 


Tube from P to Q via S, open at both ends 


i 


a 


+ 7 


Tube from P to R via S, open at both ends 


j 




/? + 7 


Tube from Q to R via S, open at both ends 


k 


a 


+ P + 1 


Forked tube joining all of P, Q, R via S, open at all three ends 


1 


a 




Surface b with large "balloon" disc attached inside the pillow 


m 


a 




Surface b with punctured torus attached inside the pillow 


n 


a 




Surface e with punctured torus attached inside the pillow 


o 


a 


+ P 


Surface h with punctured torus attached inside the pillow 


P 


a 


+ 7 


Surface i with punctured torus attached inside the pillow 


q 


a 


+ P + 1 


Surface k with punctured torus attached inside the pillow 



Table 1: The 17 vertex normal surfaces within a 4-block 



Now consider the full set of 4-blocks Bi,...,Bk', let aj,bj,...,qj denote the corresponding 
surfaces in Bi, and let aj,/3j,7j denote the corresponding boundary curves. Any normal surface in 
Xjs is a union of normal surfaces in B\, . . . , B^, and hence can be expressed as 

(Ai,i ai + . . . + Ai,i 7 qi) + . . . + (A fcj i a fe + . . . + A M7 q fe ) 

for some family of constants Ai ; i, . . . , Xk,n > 0. In this form, it can be shown 5 that the matching 
equations for reduce to the following statement: 

There is some non-negative /i£t such that, for every i, the sum A^i + . . . + \i,n q« 
has boundary ^ai + (ify + /ryj. 

In other words, the portion of the overall surface within each 4-block B{ must have boundary 
[MXi + [iffi + fiji, where [i is independent of i. 

Return now to a single 4-block with admissible vertices a, . . . , q, and let Ai a + . . . + A17 q be 
some point in the projective solution space for this 4-block. We can ensure that the corresponding 
surface has boundary of the form fia + + /ry by imposing the following constraints: 6 

A2 + A5 + A§ + A9 + An + A12 + A13 + A14 + A15 + Ai6 + A17 

= A3 + A6 + A§ + A10 + An + A15 + A17 

= A4 + A7 + Ag + A10 + An + Ai6 + A17 

This has the effect of intersecting the original projective solution space for the 4-block with two 
new hyperplanes. A standard application of the filtered double description method [9] shows that 
the resulting polytope has 18 admissible vertices, described by the following 18 normal surfaces: 
the original a with no boundary, and 17 new surfaces 7 all with boundary a + (3 + 7. Within each 
block Bi, we label these 17 new surfaces Vj 1, . . . , Vi^n- 

'The argument uses the facts that curves aj,/3»,7i surround vertices Pi,Qi,Ri respectively, and that all of these 
vertices are identified together in the overall triangulation Xk- 

6 Each line in these constraints corresponds to a section of the Boundary column in Table 1. 

7 These are the six surfaces (c + g, d + f , or j) + (b or 1), the five surfaces c + d + (b, e, 1, m, or n), and the 
six surfaces c + i, c + p, d + h, d + o, k and q. 
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Given the formulation of the matching equations above, it follows that the normal surfaces in 
Xk are described completely by the linear combinations 



Pi,i vi,i + . . . + p k ,i7 v fe ,i7 + r/iai + . . . + n k a k , 



where each Pij,r]i > and where ^ ■ pij = ^) ■ p2,j = ■ ■ ■ = Pk,j- The full projective solution 
space for X k therefore has 17 fe + k admissible vertices, corresponding to the k surfaces ai, . . . ,ajt 
and the I7 k combinations v ltjl + v 2j2 + . . . + \ k ,j k for ji,j 2 , ■ ■ ■ , jk € {1, . . . , 17}. □ 

The pathological triangulations X 1 ,X 2 ,... cover all sizes of the form n = 4k. We can generalise 
this construction to include n = 4k + 1, 4k + 2 and 4/c + 3 by replacing one of our 4-blocks with a 
single "exceptional" block. The general constructions and analyses are detailed in the full version 
of this paper, and the results are summarised in the following theorem. 

Theorem 8. For every positive 1,2, 3, 5, there exists a closed 3-manifold triangulation of size 
n whose admissible vertex count is as follows: 

n = 4k (k > 1) =► a = 17 k + k 

n = 4k + l (k > 2) => a = 581 ■ 17 fc ~ 2 + jfe + 1 

n = 4/c + 2 (jfe > 1) a = 69 • 17*- 1 + jfe ^ 

n = 4k + 3 (A; > 1) =^ a = 141 • 17 fe - x + fe + 2 

Lemma 7 proves this result for the first case n = 4k. For an extra measure of verification, 
equation (2) has been confirmed numerically for all n < 14 by building the relevant triangulations 
and using Regina to enumerate all vertex normal surfaces. 

The main result of this section is the following limit on any upper bound for a, which follows 
immediately from Theorem 8. Moreover, as we discover in the following section, there is reason to 
believe that this may in fact give the tightest possible asymptotic bound. 

Corollary 9. Any upper bound for the admissible vertex count a must grow at a rate of at least 
n(17 n / 4 ) ~ O(2.03 n ). 



5 Practical Growth 

We turn now to a comprehensive study of the admissible vertex count a for real 3-manifold trian- 
gulations. The basis of this study is a complete census of all closed 3-manifold triangulations of 
size n < 9. This is a significant undertaking, and such a census has never been compiled before; 
the paper [5] details some of the sophisticated algorithms involved. 

The result is a collection of 149 676 922 triangulations, each counted once up to isomorphism (a 
relabelling of tetrahedra and their vertices). It is worth noting that within this large collection of 
triangulations there is a much smaller number of distinct 3-manifolds, as indicated by the 3-manifold 
census data of Martelli and Petronio [26] and the author [5]. 

For each of these ~ 150 million triangulations we enumerate all vertex normal surfaces using 
the algorithms described in [7, 9]. The resulting admissible vertex counts a are summarised in 
Table 2. All computations were performed using the software package Regina [3, 4]. 

The figures that we see are remarkably small. For n = 9 tetrahedra, although Theorem 1 places 
the theoretical bound at ~ 0(29 n ), we have just 584 vertex normal surfaces in the worst case. The 
mean admissible vertex count for n = 9 is much smaller again, evaluated at just 78.49. The full 
distribution of all admissible vertex counts for n = 9 is shown in the left-hand graph of Figure 8. 
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Number of 


Number of 


Admissible vertex count (<r) 


tetrahedra (n) 


triangulations 


Mean 


Std dev 


Min 


Max 


1 


4 


2.00 


0.71 


1 


3 


2 


17 


3.94 


1.39 


2 


7 


3 


81 


5.49 


1.97 


2 


11 


4 


577 


8.80 


3.38 


2 


18 


5 


5184 


13.34 


5.49 


4 


36 


6 


57 753 


20.76 


9.21 


4 


70 


7 


722 765 


32.17 


15.29 


4 


144 


8 


9 787509 


50.20 


25.52 


4 


291 


9 


139103 032 


78.49 


42.51 


4 


584 



Table 2: Summary of admissible vertex counts for all triangulations (n < 9) 



Admissible vertex counts for n = 9 tetrahedra Growth of admissible vertex counts 




100 200 300 400 500 600 1 2 3 4 5 6 7 8 9 



Number of vertex normal surfaces Number of tetrahedra 

Figure 8: Aggregate results for admissible vertex counts 



Indeed, our pathological triangulations Xi,X2 are the worst cases for n = 4, 8 respectively, 
giving the maximum observed values ofcr=17 1 + l = 18 and <r = 17 2 + 2 = 291. More generally, 
the pathological triangulations of Theorem 8 give the maximum cases in our census wherever they 
are defined (i.e., 1,2, 3, 5). This leads us to the following general conjecture: 

Conjecture 1. For every positive n / 1,2,3,5, equation (2) gives a tight upper bound on the 
admissible vertex count a. As a consequence, we have a G 0(17 n / 4 ). 

The growth rate of a for n = 1, . . . , 9 is illustrated in the right-hand graph of Figure 8 (note 
that the vertical axis is plotted on a log scale). The growth rate of the maximum a is roughly 
17n/4 ^ 2 .03 n as suggested above; the growth rate of the average a is in the range 1.5 n to 1.6 n . 
This is just below the Fibonacci growth rate of (f> n ~ 1.62™. Indeed, if we let a n denote the 
mean admissible vertex count amongst all triangulations of size n, we find that o~ n <C o~n—i ~t~ o~n—2 
throughout our census. This leads us to our next general conjecture: 

Conjecture 2. For every n > 3, the mean admissible vertex count a n satisfies the relation a n < 
&n-i + &n-2- As a consequence, a n is bounded above by 0((j) n ) where (f> = (1 + y/b)/2. 

In particular, our census analysis gives us the following computational result: 

Theorem 10. Conjectures 1 and 2 are true for n < 9. 
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6 Conclusions 



We have pushed the theoretical bounds on the admissible vertex count a from both directions, and 
we have shown through an exhaustive study of ~ 150 million triangulations that a is surprisingly- 
small in practice. We close with a brief discussion of the implications of this study. 

Most importantly, it suggests that topological algorithms that employ normal surfaces might 
not be as infeasible as theory suggests. Hints of this have already been seen with the quadrilateral- 
to-standard conversion algorithm for normal surfaces [7], which (against theoretical expectations) 
appears to have a running time polynomial in its output size. 

In many fields, a census for size n < 9 might not seem large enough for drawing conclusions 
and conjectures. However, there is evidence elsewhere to suggest that 3-manifold triangulations 
are flexible enough for important patterns to establish themselves for very low n. For example, 
the papers [6, 29] discuss several combinatorial patterns for n < 6; these patterns have later been 
found to generalise well for larger n [5, 26], and some are now proven in general [18, 20]. 

Finally, it is clear from this practical study that the theoretical bounds on a still have much room 
for improvement. One possible direction is to incorporate the quadrilateral constraints directly into 
McMullen's theorem. This is difficult because the quadrilateral constraints break convexity, but 
the outcome may be significantly closer to the 0(17 n / 4 ) that we see in practice. 
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